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Abstract 

We present a full identification of lattice model properties with their field theoretical 
counter parts in the continuum limit for a supersymmetric model for itinerant spinless 
fermions on a one dimensional chain. The continuum limit of this model is described by 
an A/" = (2,2) superconformal field theory (SCFT) with central charge c = 1. We identify 
states and operators in the lattice model with fields in the SCFT and we relate boundary 
conditions on the lattice to sectors in the field theory. We use the dictionary we develop 
in this paper, to give a pedagogical explanation of a powerful tool to study supersymmetric 
models based on spectral flow [1]. Finally, we employ the developed machinery to explain 
numerically observed properties of the particle density on the open chain presented in [2]. 
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1 Introduction 

In the past decade a model that incorporates supersymmetry in lattice models of strongly 
interacting spinless fermions has been introduced [3] and explored for a variety of lattices in one 
and higher dimensions I31[5l[6l[71[2l[8l[9l[Il[l0l[IIl[l2l[I3] (for reviews see [H [E]). These 
models show various interesting features that are closely linked to the supersymmetry. First, 
supersymmetry provides a rich mathematical structure that allows for a considerable degree of 
analytic control in the regime where more standard perturbative techniques would fail. Second, 
it induces a delicate balance between kinetic and potential terms, which results in a strong form 
of quantum charge frustration. This so-called 'superfrustration', which typically occurs for two 
(or higher) dimensional lattices, is characterized by an extensive ground state entropy [5l [7]. 
For the one dimensional chain, the model was found to be integrable and quantum critical. In 
recent work [l2l[13], it was shown that the one dimensional model also enjoys special features, 
such as scale-free properties, which are again directly related to supersymmetry. 

In this paper, we focus on the supersymmetric model on the one dimensional chain. A 
Bethe Ansatz soltution for this model was presented in |3l [J . The continuum limit is described 
by an = (2,2) superconformal field theory (SCFT) with central charge c = 1. Here, we 
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present a full identification of lattice model properties with their field theoretical counter parts 
in the continuum limit. In particular, we identify states and operators in the lattice model 
with fields in the SCFT and we relate boundary conditions on the lattice to sectors in the field 
theory. We will see that this model forms a textbook example of how a superconformal field 
theory can be identified studying the finite size lattice model properties. Apart from finite 
size scaling of the spectrum, we will study a boundary twist, entanglement entropy and one- 
point functions. The boundary twist in the lattice model is related to a spectral flow in the 
continuum theory. The power of this type of analysis was demonstrated in [Ij where it was 
used to identify quantum criticality in various ladder models. The discussion here serves as 
a pedagogical introduction to this technique. Finally, the dictionary developed in this paper 
allows one to study various properties of the model. As an example we compute the density of 
the chain with open boundary conditions. This property was studied in [2], where its scaling 
dimension was identified numerically. Furthermore, a remarkable Z3 substructure was observed. 
The result we obtain here, both confirms and explains these observations. 

The paper is organised as follows. We first introduce the model and discuss some of the 
basics of supersymmetry. We then discuss the N = (2,2) superconformal field theory with 
central charge c = 1. In section |4j we present the full identification of this theory with the lattice 
model with closed boundary conditions. The next section provides a detailed explanation of the 
spectral flow analysis. Section [6] briefly presents some results for the entanglement entropy. In 
section [7j we extend the dictionary to the case of open boundary conditions, which is put to 
work in section [8j where the site-dependend particle density is computed on the field theory 
side. 

2 The model 

In quantum mechanics, supersymmetric theories are characterized by a positive definite energy 
spectrum and a twofold degeneracy of each non-zero energy level. The two states with the same 
energy are called superpartners and are related by the nilpotent supercharge operator. Let us 
consider an A/" = 2 supersymmetric theory, defined by two nilpotent supercharges Q and [IS] , 

g2 = (Qt)2 ^ 

and the Hamiltonian given by 

H = {Q\Q]. 

From this definition it follows directly that H is positive definite: 

= IQIV')I' + IQV)I'>o. 

Furthermore, both Q and Q'^ commute with the Hamiltonian, which gives rise to the twofold 
degeneracy in the energy spectrum. In other words, all eigenstates with an energy Eg > {) form 
doublet representations of the supersymmetry algebra. A doublet consists of two states |s), Q|s), 
such that = 0. Finally, all states with zero energy must be singlets: Q\g) = Q'^\g) = 

and conversely, all singlets must be zero energy states [16]. In addition to supersymmetry our 
models also have a fermion-number symmetry generated by the operator F with 

[F,Qt] = _Qt and [F,Q]=Q. 
Consequently, F commutes with the Hamiltonian. 
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We now make things concrete and define a supersymmetric model for spin-less fermions 
on a lattice, following [3]. The operator that creates a fermion on site i is written as c| with 
{c|,Cj} = Sij. To obtain a non-trivial Hamiltonian, we dress the fermion with a projection 
operator: P<i> = Yij next to ii^ ~ ^j^j)^ which requires all sites adjacent to site i to be empty. 
With Q = ^c|P<i> and = ^CjP<j>, the Hamiltonian of these hard-core fermions reads 

i j next to i i 

The first term is just a nearest neighbor hopping term for hard-core fermions, the second term 
contains a next-nearest neighbor repulsion, a chemical potential and a constant. The details of 
the latter terms will depend on the lattice we choose. 

In this paper we consider the supersymmetric model on a one dimensional chain. For a chain 
of length L, the Hamiltonian takes on the explicit form: 

L L 

H = Y^ [^*-i(4ci+i + 4^iCi)Pi+2\ + E(nini+2) + L-2F. 

i=l i=l 

Here Pj = 1 — n,, rii = clci is the usual number operator and F = is the total number of 

fermions. 



3 Continuum theory 

The supersymmetric model on the chain can be solved exactly through a Bethe Ansatz [3j. In 
the continuum limit one can derive the thermodynamic Bethe Ansatz equations. The model has 
the same thermodynamic equations as the XXZ Heisenberg spin chain at a specific value of the 
anisotropy parameter A. There is indeed a mapping between the supersymmetric model on the 
chain and the Heisenberg spin chain with special boundary conditions 0]. The hamiltonian of 
the XXZ chain is defined in terms of the usual Pauli matrices as 

Hxxz = \ E [«+i + ^r^r+i - ^f+i] . (1) 

i=l 

The continuum limit of the XXZ chain is described by the massless Thirring model ^17j, or 
equivalently a free massless boson <I> with action |18j 

S = ^jdxdt [{dt^f-{d,^f]. (2) 

The coupling constant g is related to the anisotropy parameter A in the XXZ chain. On the 
conformal field theory side, the coupling constant g is related to the compactification radius 
R of the free boson theory via g = 2/B?. The free boson theory is characterized by a central 
charge c = 1 and a set of highest weight states which depend on the compactification radius. 
At a compactification radius R = \/3, the conformal algebra is enhanced to an A/" = (2, 2) 
superconformal algebra [19', flH]. The (2,2) means that both the holomorphic as well as the 
anti-holomorphic fields satisfy an = 2 superconformal algebra. More specifically, the free 
boson at compactification radius R = \/3 is the simplest field theory with N = (2,2), namely 
the first in the series of minimal supersymmetric models. It turns out that the compactification 
radius R = \/3 corresponds to an anisotropy parameter of A = —1/2 (see for example j2U] ) 
which is precisely the value one obtains upon mapping the supersymmetric model onto the XXZ 
chain [5]. 

The fact that the low-energy spectrum of the supersymmetric model on the chain is described 
by a superconformal theory in the continuum limit, tells us that the model is quantum critical. 
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3.1 Superconformal algebra at R = \/3 

In an = 2 superconformal field theory |2H [22] there are three generators besides the stress- 
energy tensor: two supercharges, G~^{z) and G~{z), with conformal dimension 3/2 and a U{1) 
current, J{z), with conformal dimension 1. The N = 2 superconformal algebra is then given by 
the Virasoro algebra 

1 „ 

[Lm, Ln] = (m - n)Lm+n + Y^c(m - m)5m+nfl- (3) 
together with a C/(l) Kac- Moody algebra for the current 

\JmiJri\ — g'^^m+rijO [Lmj Jn\ — '^Jm+m (4) 

and the algebra of the supercharges 

[Lm.Gf] = (^m-r)G^+^, (5) 

[JmjG^] = ±G^_,_^, (6) 

{Gf,G^} = 2Lr+s±{r - s) Jr+s + \c{r^ -\)5r+s,o- (7) 

Here r runs over all values in Z + a, with a a real number which determines the branch cut 
properties of G^{z). For a = the theory is said to be in the Ramond sector and for q = 1/2 
it is said to be in the Neveu-Schwarz sector. 

We will now identify the supercharges and the U{\) current in the free boson spectrum at 
compactification radius R = \/3- The free boson field ^ can be decoupled into left and right 
movers: ^ = + <I>/j and the dual is defined as <I> = g{^L — ^r)- The left moving field obeys 
the following OPEs 

^l{z)^l{w) ln(2; - w), d^L{z)d^L{w) ^- — ^— ^, 

2g 2g [z - wy 

and similarly for the right moving field 

^r{z)^r{w) ~ ln(z - u7), d^R{z)d^R{w) ~ ^ 

2g 2g [z - wy 

The operators 

Vm,n = ■ exp(im<I) + m^) : , 

are called vertex operators. Here the semicolons imply normal ordering, in the following we 
will drop this notation and tacitly assume that normal ordering is taken care of. The vertex 
operators are primary fields with conformal dimensions: 

hL,R = {m±gnf/{Ag), (8) 

with m £ Z and n € Z/2. Note that we label holomorphic and anti-holomorphic dimensions 
with L and R, for left and right movers, respectively. From ^ we find that for R = \/3, and 
thus g = 2/3, the operators l±i _|_3/2 have conformal dimensions {hL,hfi) = (3/2,0) and the 
operators y±i^zps/2 have conformal dimensions {hi^hR) = (0,3/2). These four operators are the 
two left mover supercharges and two right mover supercharges. Finally, the U{1) currents of 
dimensions (1,0) and (0, 1) are proportional to d^i and d^R respectively. The proportionality 
factor follows from comparing the OPEs of d^L,R with the OPE of the U(l) current 

J{z)J{w) ~ Jl,r{z) = ±V25c/39$L,ij = ±i2/3d^L,R. (9) 

[z w) 

So these operators, together with the stress-energy tensor form an A/" = (2, 2) superconformal 
algebra. That is, both the left-moving and right-moving operators generate an M = 2 super- 
conformal algebra. 
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3.2 Spectrum 

The hamiltonian, i.e. the energy operator, is the generator of translations in the time direction. 
On the cyhnder it reads 

H = Lla + Lr^o - c/12. (10) 

It follows that for c = 1 the energy of a state is given hy E = hji + — 1/12. Just as for the 
lattice model, we find that, if we define the two supercharges 

1 1 

we can write the hamiltonian on the cylinder as 

H = Ll,o + Lr,o-^ = {G,G^}. (12) 

Note that these supercharges are defined in the Ramond sector. Indeed it turns out that the 
supersymmetric structure as we discussed it in section [2] is only fully realized in the Ramond sec- 
tor of the superconformal algebra. From the commutator of the supercharges with the Virasoro 
generators one easily verifies that the hamiltonian commutes with G and G^ . Further- 
more, from the commutation relation with the U(l) current, we find that the hamiltonian also 
commutes with F = Jl q — Jr^. The supercharges satisfy 

[F,G] = -G \f,G^] =G^. (13) 

These are precisely the relations we found for the supersymmetric model. It thus follows that 
the spectrum of the hamiltonian in the Ramond sector will be positive definite and decomposes 
into zero-energy singlets and positive energy doublets. 

3.3 Supercharges and the U{1) current 
To see the action of the supercharges 

^^ = ^^1,^3/2 and G± = y±i,^3/2. (14) 

we consider the action of the supercharge on a state Vm,n\0)- The mode expansion for G^ 
is given by 



where Z £ Z in the Ramond sector and / G (Z + 1/2) in the NS sector. Consequently, we have 

/dz 



'^'^^/^™+"+'+^/V„+i,„+3/2|0), (15) 



27rz 

where, in the second line, we used the OPE for Vi^'^i2{z)Vm,n{w^w): 

Vi,^/2{z)V^^n{w,w) ~ (z-«;)3/2™+"v;„+,_„+3/2 («;,«;). 
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Similarly, we obtain 

For the contour integrals to be well-defined the power of z or z has to be integer. Since I is 
integer (half- integer) in the Ramond (NS) sector, we find 3/2m -|- n is half-integer (integer) in 
the Ramond (NS) sector. This condition can be reformulated by saying that (_i)'"+2'^ i 
the Neveu-Schwarz sector and —1 in the Ramond sector. 

Furthermore, we see that the first mode ^ of that gives a non-zero contribution must 
obey 3/2m + n + I + 1/2 < —1, i.e. / < —3/2 — 3/2m — n. Equivalently, we find for G^f- 
I < -3/2 + 3/2m + n and for G^ ^: / < -3/2 =F 3/2m ± n. 

Let us now determine the U{1) charges of the vertex operators. The left moving U{1) charge, 
qi, is defined by the OPE of the U{1) current with a primary field ^ 

[z — W) 

and similarly for qr. Using the OPE 

z(3m + 2n) Vm,n 



d^L{z)Vm,n{w,w) 



{z — w) 



and a similar expression for the right movers, we find the U{1) charges corresponding to the 
vertex operator Vm,n to be 

qL,R = n/3±m/2. (16) 
3.4 Highest weight states and descendants 



Let us identify the highest weight states. Note that the supercharges ( 14) change m by ±1, while 
leaving (_i)'"+2n unchanged. Furthermore, by acting on a state with certain combinations of 
the supercharges, one can raise or lower n by multiples of 3 while keeping m fixed. From this, 
it follows that we need only three highest weight states per sector, since all other states can 
be generated from these states with the supercharges. For example, one can easily check that 
^o,-5/2|0) = G+ _,G2_,Vo^y,\0) and y-i,-i|0) = G^,_il^o,i/2|0>. 

In the Ramond sector, we choose the highest weight states Vo,±i/2|0) and V"o^3/2|0), since these 
are the states with lowest energy. Their conformal dimensions are /il,_r = 1/24 and /il,_r = 3/8, 
respectively, and thus their energies are E = and E = 2/3. Clearly, the same reasoning applies 
in the Neveu-Schwarz sector and the highest weight states are Vo,o|0) and Vb,±i|0) with hL,R = 
and hL,R = 1/6 respectively. Note that these are precisely the conformal dimensions of the 



first minimal model in the supersymmetric minimal series \23\ I24|. Using (16), one also readily 
verifies the agreement with the first minimal model for the U{1) charges of the highest weight 
states. 

The spectrum is generated from the highest weight states by acting with the Virasoro gen- 
erators and the supercharges. The action of the Virasoro generators on a highest weight state 
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I 1 ^r) is well known 



LL/R,o\hL,hR) 
I'L/R,n\hL, ^r) 
LL-n\hL, hR) 
LR-n\hL, hR) 



hL,R\hL, hR), 
0, 

\hL + n,hR), 
\hL, hR + n), 



with n > 0. The vacuum is defined as the state with L^^^jiolO) = 0. 

Supersymmetry implies that the zero energy states in the Ramond sector do not have su- 
perpartners. Consequently, they must be annihilated by the zero modes of the super char ges, 

we 



From the inequalities relating the mode I to m and n given in section 



3.3 



that IS G^^j^ Q 

find that indeed G^y^ qV^ -1-1/2 10) = 0. The third highest weight state in the Ramond sector, 
hL,R = 3/8, has non-zero energy, so we would expect this state to have a superpartner. In 
fact, since the left and right movers completely decouple in the continuum limit, the continuum 
theory has two N = 2 supersymmetries. Consequently, the third highest weight state forms a 
quadruplet instead of a doublet. We find that there are four states with hL,R = 3/8 and energy 
E = 3/8 + 3/8 — 1/12 = 2/3, which are all related via the supercharges 

GZ,oWlO) = ^-i,o|0), 

GZ,oG^,oW|0) = G+oG^_oT/o,3/2|0) = T^o,-3/2|0). 
A pictorial summary of the above can be found in figure [T] 
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Figure 1: The table shows hL + hR = S/im^ + 1/3^^, where (-l)"^+2" ig +1 in the NS and -1 
in the R sector. The big arrows point at the superpartners of the state with h^ + hR = 3/4. The 
symbols can be found in the 3D plot, where the energy E = hi + hR — c/ 12 is plotted against 
n and m. 



4 Relation to the lattice model 

4.1 Relation between sectors and boundary conditions 

As was noted before, the supersymmetric structure as we discussed it in section [2] is only fully 
realized in the Ramond sector of the superconformal algebra. It turns out, however, that the NS 
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is also realized in the lattice model, namely when we impose anti-periodic bomidary conditions. 
In particular, we find from numerics that the lowest energy state for the supersymmetric model 
on a periodic chain with anti-periodic boundary conditions and length L = mod 3, has a 



negative energy. A scaling analysis shows (see section 4.4) that the energy of this state is 
E = —1/12, corresponding to the NS vacuum |0). 

In the previous sections we have identified the Ramond sector of the field theory to correspond 
to the lattice model on a chain with length L = mod 3 with periodic boundary conditions. 
For the chain of length L = 3j, the fermion number in the ground state sector is fcs = L/3 = j. 
For a chain of length L = 3j±l, we would correspondingly find fas = L/3 = jibl/3. We know, 
however, that for a chain of length L = 3j it 1 the ground state has fermion number fcs = j- 
It follows that, compared to the chain of length L = 3j, the chain of length L = 3j it 1 has a 
slightly lower /higher charge density in the ground state sector. Now remember that m gives the 
charge compared to the charge in the ground state sector of the chain of length L = 3j, that is 
L/3. Suppose that instead of m = / — fas, we now write m = f — L/3. For the chain of length 
L = 3j, the two definition are completely equivalent. However, for a chain of length L = 3j it 1 
we now find that m = fas ~ L/3 = =Fl/3 in the ground state sector. The supersymmetry in 
the lattice model tells us that the chain with periodic boundary conditions corresponds to the 
Ramond sector in the field theory. Now that m takes values in Z =p 1/3, it follows that in the 
Ramond sector, which has 3m/2 -|- n E Z + 1/2, n is now integer. In fact, we find that adding 
or subtracting one site from a chain of length L = 3j corresponds in the field theory to acting 
with the operator V^i/3,_i/2- Upon acting with this operator the highest weight states in the 
Ramond sector become 

V^TiA-ilO), l^TiAolO) and VVi/3,i|0), (17) 

for the chains of length L = 3j it 1. The corresponding energies are respectively E = 1/3, and 
1/3. This agrees with the fact that these chain lengths only have one zero energy ground state. 

4.2 Lattice operators: fermion number and momentum 



From section 3.2 we know that the operator F = J^^ — Jji q satisfies the same commutation 



relations with the supercharges as the fermion number operator in the lattice model. Using 



the definition of the U{1) charges (16), we find that this operator gives Ql — Qr = fn for a 
state Vm,n\0)- If we identify the states Vq _|_i/2|0) with the two zero energy ground states of the 
supersymmetric model on the one dimensional periodic chain with length L = mod 3, which 
have fermion number /qs = L/3, we can identify the fermion number operator with 

F = Jl,o - Jr,o + fas- (18) 

In the field theory the operator that generates translations in the space direction corresponds 
to rotations on the complex plane. The momentum operator on the lattice is thus proportional 
to Lx^ Q — LfiQ. All highest weight states in the field theory have = hji, which would imply 
zero momentum. However, the two ground states of the periodic chain of length L = 3j have 
momenta po = ±7r/3 -1- irfcs mod 27r [1]. 

Using the fact that a boundary twist in the lattice model corresponds to a spectral fiow (see 
also section [s]) in the field theory, we can identify the momentum of a highest weight state with 
the sum of its U{1) charges. 

A boundary twist in the lattice model is defined by the condition that the wavefunction 
picks up a factor exp(27rza) when a fermion hops over the end of the chain, i.e. between site L 
and site 1. For a = we have periodic boundary conditions, which corresponds to the Ramond 
sector, whereas for a = 1/2 we have anti-periodic boundary conditions, corresponding to the 
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Neveu-Schwarz sector. Consequently, the boundary twist in the lattice model corresponds to a 
spectral flow in the continuum theory. 

In the following we will show, on the one hand, that the momentum of a state in the lattice 
model depends linearly on the boundary twist and, on the other hand, that the index n of a 
vertex operator Vm,n will change linearly under spectral flow. These observations will allow us 
to relate the two. 

Momentum, p mod 2tt, can be defined by writing the eigenvalues of the translation operator 
as t = e*^. The boundary twist can be implemented by replacing the term that hops a particle 
over the boundary c^i^ci+ h.c. by e'^'^*"(c|^ci+ h.c). The eigenvalues of the translation operator 
for general a then follow from 

T^\^) = e*«^e2"*"^|V') = e*P"^|V>, (19) 

where po is the momentum of {il^) for a = 0, L is the length of the system and / is the total 
number of particles in the state It follows that momentum indeed depends linearly on the 
boundary twist: Pa = Po + 2Traf/L mod 27r. 

It the continuum theory, the spectral flow is induced by the operator Vo,^. Note that it con- 
serves the fermion number, but changes the sector of the theory: (|_i)™+2n _^ (^_i'^m+2n^_'^-j2a _ 

If we now combine the fact that Pa = Po ~^ '^'^ctfcs/ L + 2Tiam/L mod 2ti for a state in the 
lattice model and n = uq + a for an operator in the field theory, we find that p is proportional 
to n. By eliminating a and using the known results for the zero energy ground states in the 
lattice model we obtain 



p = 2TTn/3 + 2'irnm/ L + fcsT^ mod 27r. (20) 

Note that the middle term is precisely 27r(/i/j — hL)/L. 

For the Ramond vacua Vo_-i-i/2|0) we can easily check this relation. We know that the two 
ground states of the periodic chain of length L = 3j have momenta po = ±7r/3 + it fas mod 2tt. 
Since the Ramond vacua have n = ±1/2 and m = 0, we find that this indeed nicely agrees with 
the equation above. Furthermore, one readily verifies that the ground state of chains of length 
L = 3j zb 1, which we identified in the previous section with the field V^i/s^O) ^-Iso has the correct 
momentum; p = + jn mod 27r. 



Finally, using the definition of the ^7(1) charges (16) we can write momentum as 

P= {ql + Qr)-^ + 2TT{hR - hL)/L + /gstt mod 27r (21) 
and we find that the momentum operator on the lattice can be expressed as 



P = {Jl,o + Jufi)-^ + {Ll,o - LRfl)2'K/L + fcsT^ mod 27r. 



(22) 



4.3 Finite-size scaling: the Fermi velocity 

The finite-size scaling of the energy depends on the boundary conditions. For (anti-) periodic 
boundary conditions, which corresponds to the cylinder on the field theory side, the scaling is 
given by [23 ES! 

Enum = 2ttEscftvf/N + 0{1/N^), (23) 

where N is the length of the finite system and vp the Fermi velocity. It follows that by comparing 
the finite size spectra with the spectrum of the field theory one can extract the Fermi velocity. In 
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this case, however we can also obtain the Fermi velocity using the mapping of the supersymmetric 
model onto the XXZ chain. For the XXZ chain ([T]) the Fermi velocity is given by 

z;f(A) = ^sin0/0, (24) 

with cos 9 = — A. The supersymmetric model maps to the XXZ chain with A = —1/2, so we 
find 6 = tt/S and the Fermi velocity of the corresponding XXZ chain is vp = (3\/3)/2. To find 
the Fermi velocity for the supersymmetric model, we note that the length N of the XXZ chain 
is related to the length L of the supersymmetric chain via N = L — f , where / is the number 
of fermions in the supersymmetric model . In the continuum limit the low energy states have 
approximately / = L/3, so N = 2L/3. Combining all this, we find for the supersymmetric 
model that 

£^num = 27rEsCFTVF/L + 0(1/1^), (25) 

with Fermi velocity vp = 'i/'2vF,xxz = (9\/3)/4. 



4.4 Finite size spectra 

In this section we analyze the numerically obtained spectra for the supersymmetric model on the 
chain with periodic and anti-periodic boundary conditions. We consider lengths up to L = 27. 

For chains of length L = 3j with anti-periodic be we find that there is one negative energy 
state with fermion number / = L/3. We identify this state with the vacuum of the supercon- 
formal field theory. The vacuum state has conformal dimensions = = 0, so its energy is 
given by -EscFT = hi + hn — c/12 = —1/12, where we used c = 1. In figure [2] we show a scaling 
analysis of the numerically obtained energies of this state for different lengths of the system. 
Remember that the scaling is given by \25\ |26] 

^num = 27rEsCFTVF/L + 0{1/L^), (26) 

where vp = 9\/3/4. Using -Escft = —1/12, we find that the energy of these states scales as 
-E'num = — 3\/37r/(8L) ~ — 2.041/L. The function that gives the best fit to the numerics is 
f{L) = a/L + 6/L2 + c/L^ with a = -2.038, b = -0.056 and c = -6.509. Clearly, the value of 
a agrees well with the theoretical value. 

Figure [2] also shows scaling analyses for other low lying levels of chains of length L = 3j. 
Clearly, the results match well with the continuum theory. Similar results are obtained for chains 
of length L = 3j±l [13]. 

In figure [3] we plot the spectrum of the periodic chain of length L = 27 as a function of 
momentum to illustrate that the analysis in the previous section allows us to identify all the 
low lying states with fields in the field theory and their descendants. To get a nice fit with the 



numerics we have rescaled the energies of the states in the field theory using ( 26 ) . Furthermore 



the momenta follow from (21). In the plot we show the vertex operators corresponding to the 
four lowest lying states explicitly. The other states are descendants of these states, which are 
obtained by acting with the supersymmetry and Virasoro generators on these states. As an 
example, table [T] shows the operators corresponding to the descendants of the highest weight 
state Vq 1/2IO). Similarly, we can identify the low energy states of chains with anti-periodic be 
and lengths L = 3j it 1 with states in different sectors of the field theory |15j . 

An interesting point is that for the periodic chain of length L = 3j, the two states that 
correspond to the fields -1-3/2 are not degenerate at finite size. The same is true for the states 
at level 1 which correspond to the fields L2,/^^_iFo,±i/2- Since the model is exactly solvable, we 
know that this must be a finite size effect and should thus vanish in the continuum limit (see 
also [271). For the two states that correspond to the fields V^o,±3/2 5 checked this explicitly 
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Table 1: In this table we show the descendants of the highest weight state Vo^i/2|0)- This is the 
zero energy state in the Ramond sector with momentum 47r/3. We order the descendants by 
their energy, also called level, in the first column. The energy is given hy E = Hl + hn — 1/12. 
In the second column we give Hl — hji, which is directly related to the change in momentum 
with respect to the highest weight state. In the third column we give all the fields at the given 
energy and momentum. Finally, in the last column we give the superpartners of the fields at 
the given energy and momentum. Between brackets we indicate the relative charge, m, of the 
superpartners. The relative charge gives the change in fermion number with respect to the 
highest weight state. 



Energy 


Momentum 


Fields 


Superpartners (charge) 








^0,1/2 




1 


1 


-^^L,-lVo,l/2 


Gi_iW = ^-W (-1) 




-1 


-i^R -1^0,1/2 


^^^,-1^0,1/2 = ^1,-1 (+1) 





9 


\'-'L-\) ''0,1/2 
-C'L,-2Vo,l/2 


Ll-iGj^ -1^,1/2 = -f'L, -1^-1,-1 (-1) 
^^,-2^0,1/2 = ^1,2 (+1) 







-^^L,-1^R ,-1^^,1/2 
^i?-l^L-l '^0,1/2 - 1^0,-5/2 


Lr-iGj^_^Vq^iI2 = _i (-1) 

-t^L -l'^5_lVo,l/2 = Ll-iVi-1 (+1) 




-2 


(-^i? -1)^^0,1/2 
-^^,-2^^0,1/2 


^i?-2^0,l/2 = y-1,2 (-1) 

Lr,-iG^_^Vo^i/2 = Li^_iVi_i (+1) 


3 


3 


(Li,_i)^y0,l/2 

^L,-l^L,-2V'o,i/2 

-^L,-3^0,l/2 


iLL,-irv-i,-, (-1) 

Ll,-2V-1,-1 (-1) 

Ll,-iVi,2 (+1) 




1 


-^^R,-l-^^L,-2Vo,i/2 
i/L.-ll^O -5/2 


(+1) 

Li,_2Vi,_i (+1) 

iR,-l^l,2 (+1) 

Lr,_iLl,_iF_i,_i (-1) 
(+2) 




-1 


-£'L,-l(-^^i?-l)^V'o,i/2 
-C'L,-l-^^R,-2Vo,i/2 


(Lk,-i)"F-i,-i (-1) 

LR,_2y-l,-l (-1) 
i^L,-l^-l,2 (-1) 

Lr,-iLl,-iVi,-i (+1) 
Gr,-2^-i,-i (-2) 




-3 


(-^^R-l)^^0,l/2 
^R-\T-'R -■2^^,112 
^R-'iVQ,\l2 


(Lr,_i)^Vi,_i (+1) 

i^R,-2-t^l,-l (+1) 
i^il,-l't"-l,2 (-1) 
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0.05 0.1 0.15 



1 /Length 

Figure 2: On the left we plot the energy versus the inverse chain length for chains of length 
L = 3j with periodic (R) and anti-periodic (NS) boundary conditions. The numerical data and 
fits are shown. The fit function is f{L) = a/L + &/L^ + c/L^. On the right, the table shows the 
results extracted from the scaling fits. 

by verifying that the energy difference between the two states goes to zero faster than one over 
the length of the system. Indeed, we find that the energy difference scales as a/L^ + b/ , with 
a = 52 and h = —86. 

Remarkably, we find that the states with the higher energy have superpartners at / = j + 1, 
whereas the states with lower energy have superpartners at / = j — 1. This is probably best 
explained via the mapping to the XXZ chain. The supersymmetric model on a chain of length 
L at 1/3 filling corresponds to the XXZ chain of length N at zero magnetization. The two chain 
lengths are related via N = L — f, where / is the number of fermions in the supersymmetric 
model. It follows that the supercharges which add or remove a fermion in the supersymmetric 
model translate into operators on the spin chain which change the length of the chain. Since the 
energy scales as one over the length, a state with more particles in the supersymmetric model, 
which has a shorter length in the XXZ chain, will thus have a higher energy. Conversely, a state 
with less particles will have a lower energy. 

5 Spectral flow 

In this section we will discuss the effect of the spectral flow operator on the states. In the 
lattice model the spectral flow operator corresponds to the boundary twist operator. This 
correspondence has proven to be very powerful in identifying critical modes in supersymmetric 
models on ladders jl]. Since the supersymmetric model on the chain is exactly solvable, we know 
that the spectral flow should correctly describe the boundary twist. We have already seen that 
the scaling behavior of the flnite size spectra nicely corresponds to the behavior one expects 
from the held theory. However, for more complicated systems extracting the scaling behavior 
can be very challenging, whereas boundary twists are easily carried out. For this reason we will 
discuss the spectral flow and boundary twist for the chain here in quite some detail. 

The spectral flow is a map between different representations of the superconformal algebra. 
The different representations are characterized by the parameter a which appears in the mode 
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Figure 3: We plot the energy versus momentum for the periodic chain of length L = 27 at 
fermion number / = 9 (filled dots). The numerically obtained energies are rescaled with the 
factor 2ttvf/L. The spectrum for the continuum theory is indicated by the open squares, 
where the number labels the degeneracy of the energy level. For the four lowest energy states 
the corresponding vertex operators are shown explicitly. In table [T] one finds the operators 
corresponding to the descendants of the highest weight state V"o^i/2|0). 



expansion of the supercharges: 



where r runs over all values in Z + a. It can be shown that the generators of the left- and 
rightmoving superconformal algebras transform as follows under spectral flow [28| 

^L/R,n - ^L/R,n + '^-Jh/R^n + g« dn,0 
•^L/R,n = Jl/R^u + 2«'^n,0 
^L/R,r ~ ^L/R,r-a 
^L/R,r ~ ^L/R,r+a- 

One can check that for a integer the algebra maps back to itself. Furthermore, for r integer and 
a = ^ the spectral flow maps the Ramond sector onto the Neveu-Schwarz sector. 

Similarly, we have the following relations for the conformal dimensions and the U{1) charges: 

h%R = h^R + aqlji + ^a^ 

1%R = + (27) 
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It follows that energy changes parabolically with a under spectral flow 



= Eo + a{qL + qR) + -a^. (28) 

If we define Q = qi — Qr and Q = qi + Qr, which are related to charge and momentum 
respectively, we find that under spectral flow 

Qa = QO; 

Qa = Qo + ya, (29) 

that is Q is invariant and Q changes linearly with a under spectral flow. 

Remember that in the lattice model we have Pa = Po + 2Traf /L mod 27r where L is the 
length of the system and / is the total number of particles. 

To compare the numerical values we obtain for the energy of finite size systems of length 
L with a boundary twist we use Snum(«) = ^TrEaVp/ L, where vp is the Fermi velocity. Using 
the linear relation between momentum in the lattice model and the twist a, we can express the 
energy as a parabolic function of momentum 

P / N OP /r , iPa-Po)QoVF , cjpq - pp)'^ V p L 
EnumiPa) = 2ttEoVf \ 

It follows that in a finite size system we should be able to fit the energies to the following curve 

EnumiPa) = a + bpa + dp^, (30) 

where the fit parameters b and d will satisfy 

_ Qqvf _ cppLvF 

f STTp ' 

d = (31) 

If we combine this with -E'num(PQ=o) = 2-kEqVp/ L, we have three equations for four parameters 
in the continuum theory: the central charge c, the energy in the Ramond sector Eq, the sum of 
the U{1) charges Q and finally the Fermi velocity vp. It follows that from the energy dependence 
on a boundary twist, one can extract c, Eq and Q as functions of the Fermi velocity. 

5.1 Spectral flow in flnite size spectra 

In this section we present the data for periodic chains of lengths up to L = 27 with a boundary 
twist. We compute the spectrum of the system for a = |, with s = 0, . . . , 7. For s = (s = 4) we 
have periodic (anti-periodic) boundary conditions. An example is shown in figure[4| for the chain 
of length L = 27 and particle number / = 9. For the low lying states one can easily see how the 
energy changes under the boundary twist. The drawn lines are parabolic fits to the energies -Enum 
as a function of the momenta pa- For the two lowest lying states, with .Enum = for Pa=o = 2it/3 
and Pa=o = 47r/3 we find the fits f{x) = a + bx + dx^ , with (o, 6, d) = (0.607,-0.435,0.069) 



and {a,b,d) = (0.579, —0.422,0.068) respectively. Using equations (31) and the usual finite size 
scaling for the energy, we find 

Eg _ -E'num iPa)L'^ /„„x 

Qa _ b + 2dpaL 

C ~ dGTTf ■ ^"^"^^ 
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Figure 4: We show the spectra (energy times chain length versus momentum) of the 27-site 
periodic chain with 9 fermions for ck = 0,|,|,|,|,|,|,|. It follows that the spectrum given 
by the black squares is precisely the spectrum plotted in figure |3] The black line connects the 
levels for different values of the twist parameter. The red line is the parabola obtained from a 
fit to the flow of one of the Ramond vacua (the state corresponding to the field ^0^-1/2) to the 
Neveu-Schwarz vacuum (the state corresponding to the field Vo,o)- The open red (blue) circles 
on this parabola correspond to a = mod 1 (a = 1/2 mod 1). The two lowest energy states 
follow the parabola very nicely as a function of the twist parameter. For the higher energy states, 



we see that there are avoided level crossings (see also the end of section 4.4) at integer values 
of a, but for intermediate values of a they still qualitatively follow the parabola. By eye one 
can also clearly distinguish the parabola's through the first descendants of the Neveu-Schwarz 
vacuum, again interrupted by occasional avoided crossings. 

For the first fit we find {Eo/c,Qo/c) = (0,-0.334) and (£^i/2/c, Oi/2/c) = (-0.083,0.000) 
and for the second fit we find {Eq/c,Qq/c) = (0,0.334) and (-Ei/2/c, (5i/2/c) = (0.254,0.675). 
It is clear that both fits correspond quite accurately with the theoretically predicted values 
of {Eq/c,Qo/c) = (0, ±1/3) in the Ramond sector and (-E'i/2/c, Qi/2/c) = (—1/12,0) and 
(-E^i/2/c, Qi/2/c) = (1/4,2/3) in the Neveu-Schwarz sector. 

Note that the two fits are almost the same, as they should since the fields V^,_i/2 and Vo^i/2 
flow into each other under spectral flow, so their energies lie on the same parabola. 

In table |3] in appendix [A] we summarize the values we extract from the parabola fits for 
(E'i/2/c, Qi/2/c) for various system sizes. It is important to note, first of all, that the values 
are quite accurate already for very small system sizes and, second of all, that we do not have 
to compare systems of different lengths. These two properties make this analysis very powerful, 
also for systems for which we do not know what the continuum limit is. 
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6 Entanglement entropy 



Entanglement entropy is often used as a measure for the entanglement between two spatially 
separated parts of the system, but it is also a powerful technique to study criticality in finite 
size systems. Let p be the density matrix of a system in a pure quantum state |^'): p = |^')(^'|. 
Let us now divide the system in two parts A and such that the Hilbert space can be written 
as Ti = Ha (^'Hb and thus l^*) = We then define the reduced density matrix of 

subsystem A as 

PA = TvBP- (34) 

The entanglement entropy is now defined as the Von Neumann entropy of the reduced density 
matrix 



Sa = -TrpAlnpA, (35) 

and equivalently for Sb- For a system in a pure quantum state we have Sa = Sb- 

For a system with a gap and thus a finite correlation length, the entanglement entropy 
typically saturates at a certain value when the size of the subsystem exceeds a certain length 
related to the correlation length. For a critical system, which can be described by a conformal 
field theory in the continuum limit, the entanglement entropy does not saturate. In fact, it 
shows a universal scaling law at a conformal critical point \29\ [30l [3T] 

S{lA) = ^HlA) + b, (36) 

where c is the central charge of the conformal field theory and 6 is a non-universal constant. 
Finally, for a one dimensional quantum critical system of finite size L, the entanglement entropy 
scales as [311 [32] 

5(/^) = ^ln(-sin(^^)) + 6, (37) 

which reduces to the expression above for L — )• oo. 

It is now clear, that if one can compute the entanglement entropy of a one dimensional 
system, this can be a very powerful way of studying the system. It can be used to, first of all, 
determine whether the system is critical or not and, second of all, if it is critical, to determine 
the central charge of the continuum theory. An important boundary condition for this method 
to work, is that one is able to determine the entanglement entropy for a subsystem larger than 
the correlation length if the system is gapped. For exact diagonalization this is clearly not 
always the case. However, if the system is studied using density matrix renormalization group 
(DMRG) methods, this condition is often met. Moreover, the entanglement entropy comes out 
essentially for free if one determines the ground state of the system with DMRG. 

For the supersymmetric model on the chain, the entanglement entropy has been studied using 



DMRG methods ||33j. The results were fitted very well by (37) and always in good agreement 
with central charge c = 1 (with errors <^ 1%) for the ground state of chains with periodic 
boundary conditions and L = 3j it 1 or with anti-periodic boundary conditions and L = 3j. The 
ground state for the other cases, that is periodic boundary conditions and L = 3j or anti-periodic 
boundary conditions and L = 3j it 1, is degenerate. To determine the central charge reliably 
in these cases one would have to construct translational invariant ground states. Finally, for 
the chain with open boundary conditions it is difficult to obtain a precise determination of the 
central charge since the entanglement entropy is plagued by very strong oscillations of period 
3. These oscillations are due to subleading finite size corrections to the entanglement entropy 
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\34:\ [35] . Similar oscillations due to subleading corrections are observed in the fermion number 
densities (see section [s]). 

Using exact diagonalization, we determined the entanglement entropy for the ground state 
of chains with periodic boundary conditions and L = 3j ± 1, for L up to 23. The systems are 
too small to get a good determination of the central charge. They are, however, in reasonable 
agreement with c = 1. Excluding the values for < 3 and Ia > L — 3, we obtain c ~ 1.05 for 
L = 22 and c 1.04 for L = 23 (see figure [s]). 

Appart from a more thorough analysis of the entanglement entropy of this system, it would 
be interesting to look at the entanglement entropy of two disjoint intervals. It has recently 
been shown that such an analysis can in principle reveal all scaling dimensions of the continuum 
theory |36] 
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(a) L = 22 and / = 7. 
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Figure 5: Entanglement entropy versus subsystems size Ia for the ground state of the chain of 
length L = 22 with / = 7 fermions (left) and the chain of length L = 23 with / = 8 fermions 
(right). 



7 Open boundary conditions 

In this section we will discuss the chain with open boundary conditions. We identify the contin- 
uum theory (the N = 2 superconformal field theory at c = 1) and the three sectors corresponding 
to the three possible chain lengths modulo three. 



7.1 Continuum theory 

For open boundary conditions we do not expect the left- and right-moving modes to decouple, so 
the simplest guess for the continuum field theory is the N = 2 superconformal field theory with 
central charge c = 1. By comparing the spectrum of this theory with numerical computations 
of finite size spectra, we concluded that this is indeed the correct guess. 
The states of the field theory are given by the vertex operators 

^„ = e(*™*/^), (38) 

where the \/3 comes from the compactification radius R = -v/3 and in the Ramond (NS) sector 
we have m G Z + 1/2 (m G Z). The conformal dimension hm corresponding to Vm is 

hm = rn^/Q. (39) 

For m = ±3 we find the supercharges, given by 

G± = e(±*^*), (40) 
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Table 2: The degeneracy at level n is given by the number of partitions p{n). The energy at 
level n is given by -B = — 1/24 + n. 



level 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


degeneracy 


1 


1 


2 


3 


5 


7 


11 


15 


22 


30 


42 


56 


77 



with conformal dimension h = 3/2. 

In the following we only consider the Ramond sector, since this is the sector that is realized 
in the lattice model. There are three highest weight states that we need to consider. They 
correspond to the primary fields: V_i/2, Vi/2 All other states are generated by the 

supercharges and the Virasoro algebra. The fields V_i/2 and V1/2 both have conformal dimension 
h = 1/24. Since the hamiltonian is given hy H = Lq — c/24, they correspond to zero energy 
states. The state V3/2IO), however, has energy E = 9/24—1/24 = 1/3 and also has a superpartner 

{v-3/2m. 

Since the supercharges change m by ±3, we infer that m/3 corresponds to the fermion number 
in the lattice model. It thus quickly follows that the three sectors; m = 1/2, 3/2, 5/2 mod 3, are 
related to the three sectors in the lattice model with chain lengths L = 0, 1,2 mod 3. Indeed 
we know (see e.g. jUj) that chains of length L = 0,2 mod 3 have one zero energy ground states 
and chains of length L = 1 mod 3 have all energies larger than zero. It follows that chains with 
length L = 1 mod 3 correspond to the sector with m = 3/2 mod 3. To identify the sector of the 
other two chain lengths, we look at the first excited state and its superpartner. The first excited 
state is given by L_iV-|-i/2|0) and the respective superpartners are G^iV±i/2\0) = V^:f5/2|0). One 
easily checks that the superpartners indeed have energy E = - 1/24 = 1. The difference 

is that one occurs at / = fas + 1 and the other at / = fcs — 1- If we compare this with the 
finite size spectra we quickly conclude that open chains with length L = 3j correspond to the 
sector with m = 5/2 mod 3 and chains with length L = 3j — 1 correspond to the sector with 
m = 1/2 mod 3. 

Finally, we find the following relation between fermion number /, chain length L and the 
quantum number m: 

/ = /-L/3 = (m + l/2)/3. (41) 

For m = —1/2 this relation gives / = 0, which agrees with f = j and L = 3j. For m = 1/2 this 
relation gives / = 1/3, which agrees with f = j and L = 3j — 1. Finally, for L = 3j + 1 the 
two lowest energy states are found at / = j and f = j + 1, which matches with m = —3/2 and 
m = 3/2. 

The spectrum per sector now simply follows from the highest weight states and their descen- 
dants. The character formula is given by 

Xh{q) = q^/vii) 

= q'^"'Vr,{q), (42) 

where, as usual, ri{q) = q'^^'^'^\Yk=iO- ~ q^) and q = e^'^''"'^ . From this formula we obtain the 
degeneracies at each level. For the first few levels the degeneracies are summarized in table [2} 
The energy of the n-th level is given hy E = m? /Q — 1/24 + n. 

7.2 Finite size spectra 

We now compare the continuum spectrum to the numerics. We have performed a scaling analysis 
for the numerically obtained energies for chains of lengths up to L = 23. The low lying levels 
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are nicely fitted with the function f{L) = a/L + + c/L^ (for details see [15]). For open 
boundary conditions the scaling is given by 



^num = ttEscftvf/L + 0{l/L^), (43) 

where the Fermi velocity is given by vp = 9\/3/4. It follows that the energies in the continuum 
limit can be extracted from the fits via E = a/{'KVF)- As an example we show the continuum 
limit spectrum that we extracted in this way for the chain of length L = mod 3 in figure |6(b)[ 
which should be compared to the continuum spectrum for m = 5/2 mod 3 plotted in figure 



6(a) (similar results are obtained for L = 1, 2 mod 3 and can be found in [15| ) 
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(a) The energy levels E = m^/6 — 1/24 + n are 
plotted versus m for m = 5/2 mod 3 which corre- 
sponds to chains of length L = 3j. On the horizon- 
tal axis we also indicate the shifted fermion number 
/ = f—L/3 = {m+l/2)/3. The level corresponding 
to the highest weight state and the levels that are 
generated from this field by the supercharge opera- 
tors, that is the levels with energy E = m^/6— 1/24, 
are indicated by a thick bar. The descendants are 
indicated by thinner bars. The labels indicate the 
degeneracy of the levels. 




-1 
Shifted fermion number (f-L/3) 

(b) The numerically fitted values of the energy are plot- 
ted versus the shifted fermion number / = f — L/3 
for chains of length L — 3j. The fit values are ob- 
tained by fitting the numerically obtained energies as 
a function of the chain length L with the function 
f{L) — a/L + b/L'^ + c/L^. The energy is then given 
hy E — a/{nvF)- The labels indicate the number of 
overlapping data points. 



Figure 6: On the left we show the theoretically predicted spectrum and on the right the spectrum 
obtained from fits to the numerically obtained spectra. These spectra are for chains of length 
L = 3j. 

Clearly, for the first few levels, we find a very nice agreement with the theoretically obtained 
spectra. For the higher levels, the fits are not very reliable. We indicate two reasons for this. 
First of all, there is a large degeneracy of these levels in the continuum limit. However, in the 
finite size spectra the degeneracies are not realized, since for finite size matrices the eigenvalues 
tend to spread. Second of all, since there are corrections of order with p >2, there can be 
level crossings as a function of the length. Numerically, however, we simply connect the n-th 
level at different lengths, so we do not take possible level crossings into account. The latter 
argument also explains why supersymmetry may appear to be broken in the spectra extracted 
from the fits. Clearly, this is not the case in the original numerical spectra. 

Using density renormalization group methods, one can obtain the low lying levels for much 
larger system sizes. For the 120-site chain with open boundary conditions it was found that [33] 
the spectrum and level degeneracies are, at least up to level 7, in excellent agreement with the 
continuum theory in the sector with m = 5/2 mod 3. 
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8 One-point functions 



In [2] Beccaria and De Angelis discuss the supersymmetric model on the chain with open bound- 
ary conditions and length L mod 3 = 0. In particular, they use non standard number theoretical 
methods to obtain exact expressions for the ground state wave function, on the one hand, and on 
the other hand, they study the finite size scaling behavior of some simple correlation functions 
using exact diagonalization. Their results for the one point-function (n^) = (V'o|c|.Cfc|V'o) can be 
summarized as follows. They find that (n^) has a clear Z3 substructure. The one-point functions 
(nfc,fc=i mod 3) and {nk^k=o mod 3) are not symmetric under k ^ L — k. The one-point function 
{nk^k=2 mod 3) is symmetric under this map and shows a very different behavior from the other 
two. For the different branches they extract the following finite size scaling behavior 

• (nfc) - 1/3 = /+ ((fc - k+)/L^ L-" for k mod 3 = 2, 

• {nk + nk+2) - 2/3 = /_ (^{k - k^)/L^ L-" for k mod 3 = 1, 

where L = L/3 + 1, A;± = (L ± l)/2. They obtain the best fit for v = 0.33(2). 

In the following we will use the observed Z3 substructure to propose an identification of the 
one-point functions with expectation values of operators in the superconformal field theory. 

Let us first identify the Z3 operator in the superconformal field theory on the cylinder. 
Remember that the boson is compactifiec|^ <I> = <1> + 27r and $ = $ + 27r. It follows that the 
operator T that acts as fohows T : ^> H> ^> + 27r/3 and T : $ H> l> + 27r/3 satisfies = 1. We 
now consider the action of this operator on the vertex operators Vo^-ti = e^** 

rFo,±i = e±<*+2V3) ^ g±2.V3y^^^^_ (44) 

Furthermore, we trivially have TVo,o = Vb,o and clearly Vo,o is just the identity. It follows that 
Vb,o and Vb,±i are eigenfunctions of the Z3 operator with eigenvalues ujq = 1 and oj±i = e='=27ri/3 
respectively. 

In the lattice the one-point function (n^) shows strong oscillations with period three in k. 
It follows that we can identify the inverse translation operator which sends k to k + \ as the Z3 
operator]^ We thus find that the functions, 

uJiUk-i + nfc + ^Tifc+i, (45) 

are eigenfunctions of the inverse translation operator with eigenvalues uJi = g^'^*'/^. In the 
following, we will assum^ k mod 3 = 2. 
Combining these observations we find 

^oVb,o = nk-i + nk + Uk+i 

where and A±i are constants. From the fact that the ground state has filling 1/3, we 
immediately find that ^0 = 1- Solving the above equations for the we obtain 

3nfc = ^0,0 + ^1^0,1 + A_iVo,-i 

3nfc+i = Vb,o + e^^'^'^AiVo,! + e-''^''^ A^iV^^^i. 

^The normalisation in the Vertex operators is chosen such that the compactification radius R = drops out. 

^Note that there is an ambiguity here: we could just as well have identified the translation operator, which 
sends to fc — 1 with the Z3 operator. However, this ambiguity is fixed by comparison with the numerical results 
of [2]. 

^There is again an ambiguity here, since we could also take k mod 3 = 0, 1. This choice, however, is again 
fixed by comparison with the numerical results of [2] (see also the end of this section). 



21 



In the following, we compute the expectation values of these operators on the cylinder and 
on the strip, which corresponds to closed and open boundary conditions respectively. For closed 
boundary conditions we compare our findings to analytical results presented in p^. For open 
boundary conditions, we show how the expectation values can be computed analytically by 
mapping the strip onto the plane and introducing the mirror images to ensure the boundary 
conditions are preserved. The formulae we obtain are in nice agreement with the finite size 
scaling behavior found in [2]. An obvious follow-up on this work, is to extend it to two-point 
functions (see also and chains of length L ^ mod 3. 



8.1 Open boundary conditions 

The expectation values can be computed on the plane, where the correlator of vertex operators 
reads 

{Vm^nA^UZ,) . . . Vm,,n,izk,Zk)) = H^Zi - Zjf''^-^^- 1 \zi - Zjf^^-'^n,l\ (45) 

Remember that for the vertex operators Vm,n we have ol^j? = m it 2n/3. 

To illustrate how we go from the strip to the plane we compute the expectation value of the 
vertex operator Vo,i in the Neveu-Schwarz vacuum. To compute this we have to do two steps. 
First we map the strip to the upper half plane using a conformal mapping, z = e'^*"'/^, where 
z and w correspond to coordinates on the plane and the strip, respectively. Then we introduce 
a mirror image of the system on the upper half plane in the lower half plane, thus defining the 
system on the full complex plane, where we can compute the correlator. The first step gives 

{0\Vo,i{w,w)\OUnr> = (^)'"(||)'"(0|yo,i(^,^)|0)uHP 



( X 1/6 /-TTt \ 



TTt \ 1/6 f —Tri_ _ 

"""UHP, 



where we used the fact that the conformal dimensions of Vo,i are hi = hn = A correlator 
in the upper half plane can be computed using the image-technique. Choosing the boundary 
condition 5<l>|boundary = 0, that is, there is no current flow across the boundary, the image- 
technique for a vertex operator gives 

(0|K.,n(^,:2)|0)uHP = (0|y3a^/2(^)^3a«/2(^*)|0>> 

and, in particular, 

(0|yo,n(^,^)|0)uHP = mn{z)V-n{z*m, 



where the vertex operator with one index is purely holomorphic (see (38)). In the end we set 
z* = z. We thus find 

(0|yo,i(^,^)|0)uHP = (0|Fi(z)y_i(z*)|0) 

= {z-z*)~^l^. (47) 

Combining both steps and using w = x + zi, we find 

(0|yo,i(x,t)|0)strip = {^-^^'\-zzfl\z-z)-^l^ 

/7ri\l/3 g7r*{«,-«7)/6L 

= (-l)^/«(^)''%in-V3(,,/^). 
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Now remember that we identified the ground state of the chain of length L = 3j and open 
boundary conditions with the Ramond vacuum l^_i/2|0). To compute an expectation value we 



define the in-state as | ii) = lim^ 
These definitions ensure that 

{R\R) 



V^i/2{z)\0) and the out-state as {R\ = limz^oo{0\Vi/2{z)z^^'^'^ . 



lim (0|yi/2(^i)zl/'V_i/2(^2)|0) 



lim 

22^0,21— >00 
1. 



1/12, 



Zl- Z2 



-1/12 



It thus follows that 

{R\Vo,iix)\R) strip 



lim 

22— ^0,Zl— !>00 \ L 

in 



7ri\i/3 



izz)'/^{0\V^/2{zl)zl^'^Vo,liz,z)V.y2{z2)\0) 



UHP 



lim 

22— 5>0,Zl— 5>CXD V L 

-7ri\i/3 



'")'^'(zz)V6(o|yi/2(z0^i'/'Vi(z)y_i(z*)y_i/2(^2)|o) 



(-l)Va(-)-,l/3(,_,) 



^^-1/3 



In the last step we used ( 46 ) . Equivalently, we find 



(i?|^0,-l(x)|i?)stHp = {-l)'/'{^fz'/Hz-z)-'/'. 

Combining all the above, and setting Ai = A_i = A, we obtairj^ 



\2LJ sin^/^TTx/L) 
1 + 2A (—] ^'^ sin(7r(x - L)/3L) 



\2LJ sinV3 (ttx/L) 
^ _ i—V'^ sin(^3;/3L) 
\2l) sinV3(7rx/L)' 



(48) 



These equations clearly reproduce the observed scaling behavior. Comparison with the numerics 
suggests that we should choose A 0.77. Finally, we argue that we should identify the width 
of the strip L with Lc + 3, where Lc is the length of the chain. The way we understand this, 
is that the open chain can be obtained from a periodic chain with three sites extra, by pinning 
one particle to a certain site. Due to the hard-core character of the particles, the neighboring 
two sites must be empty. One thus effectively takes out three sites from the system and is left 
with an open chain. Finally, this implies that the sites -1 and Lc + 2 are identified with the 
boundaries of the strip: x = and x = L. Consequently, an arbitrary site p of the chain should 
be identified with the point x = p+1 on the strip. In figure 7(a)| we show the one-point functions 
for a chain of length Lc = 30, this plot can be compared directly to the data presented in figure 
1 in [2]. In figure [7(b)] we plot the finite size scaling functions to be compared with figure 2 in 
[2]. The agreement with the data of Beccaria et. al. is quite convincing, however, it seems to get 
poorer upon approaching the boundary of the system. In particular, the value for (ni), which 
they obtain analytically using sophisticated number theoretical methods, cannot be reproduced 
from the field theory. 



*To obtain these results, that are in good agreement with numerical observations, we write (— 1) 



1/6 



where we use the ± sign for Vb,±i respectively. This ensures that the results for Vb,±i are related via complex 
conjugation. 
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site number (p) x/L x/L 

(a) One-point functions for chain length Lc ~ 30. (b) Finite size scaling functions versus the nor- 

malised position on the strip (which is related 
to the site number, p, via x = p + 1). 

Figure 7: On the left we show the one-point functions for a chain of length Lc = 30, this plot 
can be compared directly to the data presented in figure 1 in [2j . On the right we plot the finite 
size scaling functions, {{uk) — 1/3) and {{uk+i + nk+2) — 2/3) (L/3)^/'^, to be compared 
with figure 2 in [2]. 



We can use the scaling functions to compute Fj = X^fc fc mod 3=j('^'=)' which is the average 
fermion number in the three branches distinguished by the site number modulo three. Replacing 
the sum by an integral and using that 

cos(7r(x - l/2)/3) 
Jo sin^/^(7rx) 

sin(7rx/3) 
Jo sin^/^(7rj;) 

we find 

F2 = L/9 + A^(^)'^' |2V3«F/3 + 0.52F2/3, 

F, = Fo = L/9-A^(^)'^'|2-2/='«F/3-0.26F2/^ 

where we used A ~ 0.77 and F = L/3. Consequently, the oscillation in occupation number with 
period three in the site number, does not lead to a difference in the average fermion number in 
the three branches. However, the subleading term goes as F~^^^, which goes to zero only very 
slowly for F — ?■ 00. 

On a qualitative level the Z3 substructure can probably be interpreted as follows. Remember 
that on the periodic chain we have two ground states, which are plane waves with opposite 
momenta differing by 27r/3. It seems that for open boundary conditions these two ground states 
combine into a standing wave, which explains the observed density fluctuations. 

8.2 Closed boundary conditions 

For the chain with periodic boundary conditions and length L = mod 3 there are two ground 
states. The two ground states with momenta ib7r/3 -|- fcs'^ mod 2tt correspond to |A;i) = 
lim^^o ^o,i/2('2)|0) and \k2) = lim^_j.o ^o,-i/2(-2)|0) respectively. Using the same techniques as 
above we can compute the one-point functions using a conformal mapping from the cylinder to 
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the plane. We find 



{ki\VQ^^i\kj) 



{ki\yofi\kj) 



(A;i|^o,i|^i) 




1 for i = i 
otherwise 



otherwise 



otherwise 



{'^Y'^ fori = 2andj = l 



{"^f^ fori = landj=2 



If we define the states |ib) = (|A;i) it \ k2))/^/2, we obtain 




1/3 



(±|nfc|±) 



A 



1/3 



(±|nfc+i|±) = (±|nfc_i|±) 



A 



These expressions can be compared with the results of [13j, where they find {+\nk\+) — {—\nk\—) = 



which is in good agreement with the value {A ~ 0.77) we obtained in the previous section. 

9 Conclusions 

We have presented a full dictionary that relates the supersymmetric model for lattice fermions 
on the ID chain to the N = (2, 2) superconformal field theory with central charge c = 1. As an 
example we have shown how this thorough understanding of the continuum limit can be employed 
to compute properties such as the particle density. It is remarkable that the continuum theory 
can reproduce the site dependent oscillations in the density that are observed for finite systems. 
An obvious extension of this work is to analyze two-point functions. Furthermore, this work 
could be extended to the generalized supersymmetric models that are labelled by an index k [3] . 
The generalized supersymmetric models are found to have a quantum critical point described by 
the fc-th superconformal minimal model. A similar dictionary relating the k-th. supersymmetric 
model with the fc-th superconformal minimal model could thus be developed. Finally, the study 
of the supersymmtric model on the chain away from the critical point |12i [T3l [37] is likely to 
benefit from the detailed understanding of the model at the critical point presented here. 
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. From this we extract 
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A Results from spectral flow analysis 



Table 3: In this tabic we summarize the values we extract from the parabola fits for E112/C and 
Qi/2/c for various system sizes. Note that for each chain length we give two pairs of extracted 
values. They correspond to two different highest weight states in the Neveu-Schwarz sector. 
The values are extracted from a fit to the flow of that particular highest weight state from the 
Ramond to the NS sector. To be precise, for the chains with length L = 3j the middle two 
columns arc extracted from a fit to the flow of the state corresponding to the field Vq-\I2 to the 
state corresponding to the field Vo,05 whereas the last two columns are extracted from a fit to 
the flow corresponding to V^) 1/2 — > Vb,i. For the chains with length L = 3j it 1 the middle two 
columns follow from V^i/3,0 ^ ^i/3,i/2) and the last two columns from V^xjz^-x ^^=1/3,-1/2- 
The theoretically predicted values for L = 3j arc (£'i/2/c, Qi/2/c) = (—1/12,0) and (1/4,2/3) 
for the first and second column respectively. For L = 3 j ± 1 the theory predicts (1/12, 1/3) for 
the first and (1/12, —1/3) for the second column. 



chain 


lermion 










length 


number 


7—1 / 

E1/2/C 


Q1/2/C 


7—1 / 

E1/2/C 


Q1I2IC 


6 


2 


-0.085 


-0.004 


0.339 


0.846 


9 


3 


-0.084 


-0.002 


0.285 


0.738 


12 


4 


-0.084 


-0.001 


0.270 


0.706 







-U.Uo4 


n nm 
-U.UUi 


U.ZOo 




18 


6 


-0.084 


0.000 


0.259 


0.685 


21 


7 


-0.083 


0.000 


0.257 


0.680 


24 
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-0.083 


0.000 


0.255 


0.677 


27 
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-0.083 


0.000 


0.254 


0.675 
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2 


0.086 


0.338 


0.138 


-0.560 


8 


3 


0.084 


0.335 


0.101 


-0.407 


11 


4 


0.084 


0.334 


0.093 


-0.372 


14 
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0.084 


0.334 


0.089 


-0.357 


17 
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0.084 


0.334 


0.087 


-0.350 


20 
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0.083 


0.334 


0.086 


-0.346 


23 
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0.083 


0.334 


0.086 


-0.343 


26 
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0.083 


0.334 


0.085 


-0.341 
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0.085 


0.336 


0.108 


-0.436 


10 


3 


0.084 


0.335 


0.095 


-0.381 


13 
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0.084 


0.334 


0.09 


-0.362 


16 
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0.084 


0.334 


0.088 


-0.352 


19 
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0.084 


0.334 


0.087 


-0.347 


22 


7 


0.083 


0.334 


0.086 


-0.344 


25 


8 


0.083 


0.334 


0.085 


-0.342 
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